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Abstract. We find all the XX spin chains with perfect state transfer (PST) that are 
connected with the dual -1 Hahn polynomials Rn{x; a, /3, N). For N odd we recover 
a model that had already been identified while for N even, we obtain a new system 
exhibiting PST. 
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1. Introduction 

The transfer of a quantum state is said to be perfect if there is probability 1 of finding 
after some time at an end point the state introduced as input at an initial site. It has 
been realized that perfect state transfer (PST) can be achieved in inhomogeneous XX 
spin chains in particular, provided the 1 -excitation energy eigenvalues satisfy a simple 
spacing condition. (For reviews, see [1], [2]-) It has further been shown recently [11] 
that a unique (up to trivial rescaling) XX spin chain with nearest neighbor interactions 
and with PST, corresponds to each such spectrum. This matching proceeds through 
the association of families of orthogonal polynomials in a discrete variable to XX spin 
chains with PST. The 1-excitation eigenvalues thus intervene as the orthogonality grid 
points. We have provided an algorithm to construct the XX Hamiltonians from the 
spectral data and have shown how to obtain different models with PST from a parent 
system with that property by performing spectral surgery [TT] . 

The exact solvability of the PST models is intimately related to the characterization 
of the associated orthogonal polynomials. Among the Dunkl or -1 orthogonal 
polynomials that we have discovered lately [8], one family, namely that of the dual 
-1 Hahn polynomials, interestingly relates to PST. We shall here obtain the spin chains 
with PST to which they are associated and find in the process, a new exactly solvable 
model. 

The outline of the paper is the follows. The dual -1 Hahn polynomials and their 
relevant properties are recalled in section 1. The relation between perfect state transfer 
in XX chains and orthogonal polynomials theory is reviewed in section 2. Finally, 
depending on the parity of A^-the number of sites minus one - two spin chains with PST 
are identified and analyzed in section 3. The one for odd had already been reported 
in [5] and further studied in connection with dual Hahn polynomials in [6]. The one for 
even A^ had escaped attention so far. We conclude the paper by showing how these two 
models are related. 

2. Dual -1 Hahn polynomials 

Dual -1 Hahn polynomials Rn{x; a, /3, A^) were introduced in [9] as g = —1 limits of the 
dual g-Hahn polynomials. 

These polynomials satisfy the 3-term recurrence relation 




(2.2) 



(2.1) 
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and 

_{2N + l-a-P if n even 

\_2iV-3 + a + /3 if n odd' ^ ^ 

In compact form we have 

Un = A[n]^[N - n + l]^, hn = 2{[n]^ + [N -n]^) + l-a- (3, (2.4) 

where 

, /3-A^-l a- 

^= 2 = 2 

and 

H, = n + Ml-(-ir), (2.5) 

are the so-called "/x- numbers" which appear naturally in problems connected with the 
Dunkl operators [1]. 

It is seen that Uq = u^^i = as required for finite orthogonal polynomials. The 
positivity condition m„ > 0, n = 1, 2, . . . , iV is equivalent to the conditions 

a> N, P> N. (2.6) 

Let us define the Bannai-Ito (BI) grid 

J —a — (3 + 2s + 1 if s even, 
" \ a + /3 - 2s - 1 if s odd ■ ^ ' ' 

The polynomials Rn{x) are orthogonal on the + 1 points yo,yi, . . . ,yN of the BI grid 

N 

WsRn{ys)Rm{ys) = H0U1U2 ...Un (5„m, (2.8) 

where the discrete weights are defined as 

,, A-N/2), (l-a/2),(l-a/2-/3/2), 
^ ^ s\ (l-/3/2),(iV/2 + l-a/2-/3/2), ^"''^ 



and 



_ ^ (-iV/2).+i (1 - a/2)s{l - a/2 - /3/2), 
^"^^ " ^"'^ (l-/3/2).(iV/2 + l-a/2-/?/2).+/ ^''^^^ 

The normalization coefficient is 

y 2 Jn/2 /o ii\ 

Assume that a = N + ei, f3 = N + e2, where ei^2 are arbitrary positive parameters. 
(This parametrization corresponds to the positivity condition for the dual -1 Hahn 
polynomials.) Then it is easily verified that all the weights are positive Wg > 0, s = 
0,1,..., AT. 

Moreover, the spectral points ys are divided into two non-overlapping discrete sets 
of the real line: 

{1 - 5, -3 - (5, -7 - (5, ... , -2N + 1-5} 
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and 

{l + 5,5 + 5,9 + 5,...,2Ar-3 + 5}, 

where (5 = ei + 62 > 0. The first set corresponds to Us with even s and contains 1 + N/2 
points; the second set corresponds to ys with odd s and contains N/2 points. 
When = 1, 3, 5, ... is odd we have the BI grid Hs 

\a + [5 + 2s + l if s even , ^ ^, 

y^ = \ « o 1 -f AA 2.12 

[ —a — p — 2s — 1 if s odd 

and the recurrence coefficients 

r 4:n{N + 1 — n) if n even 
" t4(a + n)(/3 + iV+ 1 -n) if n odd ^^'^^^ 

and 

-1 — a + B if n even 

(2 14) 

-1 + a - /3 if n odd 
In compact form we have 

Un = 4:[n]^[N -n + 1]^, 6„ = 2{[n]^ + [N- n]^,) -2N-l-a-(3, (2.15) 

with C = a/2, 7] = f3/2. 

Note that the parameters ^ and rj have different expressions in (12. 4p or (12.151) . 
The positivity condition m„ > 0, n = 1, 2, . . . , iV is equivalent either to 

a > -1, (3 > -1 (2.16) 

or to a < — A^, (3 < —N. It is sufficient to consider (I2.16p . 

The polynomials Rn{x) are orthogonal on the set of + 1 points ys 

N 

^WsRn{ys)Rm{ys) = H0U1U2 ■ ■ ■ Un 5nm, (2-17) 
s=0 

where the discrete weights are defined as 

= ( .y {-{N-l)/2)s (l/2 + «/2),(l + a/2 + /3/2), 

^ ^ s\ (l/2 + /3/2),(Ar/2 + 3/2 + a/2 + /3/2)/ • ^ 



and 



-( (l/2 + a/2),+i(l + a/2 + /j/2), 

^ ^ (l/2 + /3/2),+i(Ar/2 + 3/2 + a/2 + /3/2)/^ • ^ 



The normalization coefficient is 
(1 + ^) 

° = (p±l\ ■ ^2-2°^ 

V 2 ^(Ar+l)/2 



Assume that a = — 1 + ei, (3 = — l + e2, where ei,2 are arbitrary positive 
parameters. (This parametrization corresponds to the positivity condition for the dual 
-1 Hahn polynomials for A^ odd.) Then it is easily verified that the weights are positive 
Ws>0, s = 0, 1, . . . , A^. 
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Again, the spectral points Us are divided into two non-overlapping discrete sets of 
the real line: 

{-1 - 5, -5 - 5, -9 - 5, ... , -2N + 1-5} 

and 

{-l + 5,3 + (5,7 + 5,...,2Ar-3 + 5}, 

where 5 = ei + 62 > 0. Both sets contain (A^ — l)/2 points. 

The dual -1 Hahn polynomials can be explicitly expressed in terms of the ordinary 
dual Hahn polynomials; depending on whether and n are even or odd, the formulas 
take different forms (see [9] for details). 

In [To] it was shown that the -1 dual Hahn polynomials appear as Clebsch-Gordan 
coefficients for the irreducible representations of the algebra s/_i(2) which is the q = —1 
limit of the quantum algebra slq{2). The algebra s/_i(2) is generated by 4 operators 
Jo, J+, J-, -R satisfying the commutation relations 

[Jo,J±]=±J±, [i?,Jo] = 0, {J+,J_} = 2Jo, {i?,J±} = 0, (2.21) 

where [A, B] = AB — BA and {A, B} = AB + BA. The operator R is an involution 
operator, i.e. it satisfies the property 

= I. (2.22) 

The Casimir operator Q, commuting with Jq and J± is [ID] 

Q = J+J-q-''^ - _ ^^(^ _ {q'^'-' + q''^). (2.23) 

The algebra s/_i(2) admits a nontrivial addition rule (i.e. coproduct) 

Jo = Jo O / + / (g) Jo, J± = J±(g)i? + /(g) J±, R = R^R. (2.24) 

such that the operators Jo, J±, R again satisfy relations fl2.2ip . 

The Clebsch-Gordan coefficients (CGC) of sZ_i(2) arise as overlap coefficients 
between two canonical bases in the coproduct representation space. For details on 
the identification of the CGC with -1 dual Hahn polynomials see jlOj . 



3. Perfect state transfer and orthogonal polynomials 

Consider the spin chain with Hamiltonian 

7V-1 N 
«=0 1=0 

where J; > are the constants coupling the sites / — 1 and / and hi are the strengths of 
the magnetic field at the sites / (/ = 0, 1, . . . , A^). The symbols erf, , erf stand for the 
Pauli matrices. 

Introduce the basis vectors 

\en) = (0,0,...,1,...,0), n = 0,l,2,...,Ar, 
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where the only 1 (spin up) occupies the n-th position. In that basis, the restriction J 
of H to the one-excitation subspace is given by the following (A^ + 1) x (A^ + 1) Jacobi 
matrix 

/ 6o Ji ... \ 
J Ji 6i J2 ... 



\ ... Jn bN J 
Equivalently, we have. 

ll^n+l) + bn\en) + Jn\en-l)- (3.2) 

The boundary conditions 

Jo = Jn+1 = (3.3) 

are assumed. 

Let s = 0, 1, . . . , denote the eigenvalues of the matrix J. They are all real 
and nondegenerate. They are labeled in increasing order, i.e. xq < xi < X2 < . . .xn ■ 

To the Jacobi matrix J one can associate the monic orthogonal polynomials Pn{x) 
defined by the 3-term recurrence relation 

Pn+l{x) + hnPn{x) + UnPn-l{x) = xPn{x), 

n = 0,l,...,N, P_i = 0, Po(x) = 1, (3.4) 

where m„ = > 0. 

P/v+i(x) is the characteristic polynomial 

Pn+i{x) = (x - xo){x -xi) ...{x - xn)- (3.5) 

The orthogonality relation reads 

N 

(3.6) 

s=0 

where 

hn = U1U2 ...Un 

and the discrete weights Wg > are uniquely determined by the recurrence coefficients 

Perfect state transfer (PST) occurs ^ if there exists a time T such that 

e'^'leo) = e^'^M, (3.7) 

where is a real number. This means that the initial state |eo) evolves into the state 
I Cat) (up to an inessential phase factor e^"^). 

It is known that the PST property is equivalent to the two conditions [2]: 

(i) the eigenvalues Xg satisfy 

vr 

Xs+i - Xs = j;Ms, (3.8) 
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where Ms are positive odd numbers. 

(ii) the matrix J is mirror-symmetric RJR = J, where the matrix R (reflection 
matrix) is 

/ ... 1 \ 
... 1 

R 



\ 1 ... / 

Property (ii) is equivalent to the property [H] 

(ii') the weights Wg (up to a normalization) are given by the expression 
4. XX spin chains with PST and and dual -1 Hahn OPs 

We now wish to determine the circumstances for which the Jacobi matrix J 
corresponding to the -1 dual Hahn polynomials will possess the PST property. 

Crucial is the mirror symmetry condition (ii). In terms of the recurrence coefficients 
it means 

UN-n+l = Un, bN-n = bn (4.1) 

for all n = 0, 1,...A^. 

Consider first the case of odd A^. In this case formulas fl2.15p immediately imply 
that conditions (14.11) hold iff ^ = 77, or, equivalently, iff a = /3. Under this condition we 
have 

«n = 4[n]„/2[A^-n + l]„/2, = -1. (4.2) 

The positivity condition for the measure is a < —1. 

The spectrum Xs of the Jacobi matrix J coincides with the BI grid fl2.12p and 
consists of two uniform subgrids G- and (7+ (with a step of 4 between the neighbor 
points). 

The subgrid G- consists of the (A^ + l)/2 points a;o < xi < . . . < X(Ar_i)/2: 

G_ = {-2N + 1 + 5, -2N + 5 + 6,..., -1-6}. (4.3) 

The subgrid G+ consists of the (A^ + l)/2 points X(^n+i)/2 < X(n+3)/2 < ■ ■ ■ < xn- 

G+ = {-1 + 6,4 + 6,. ..,2N -3 + 6,} (4.4) 

where 6 = 2{a + 1) > 0. There is a gap of length 26 = 4(a + 1) > between these 2 
subgrids. 

Condition (i) in this case is equivalent to the restriction 
Mo 

« = ^, M2>Mi, (4.5) 
where M2 is even and Mi is odd (M2 and Mi are assumed to be coprime). 
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The Jacobi matrix with coefficients fl4.2p and with the PST property, was first 
introduced in [5]. In [6] the corresponding orthogonal polynomials Pn{x) were related 
to the dual Hahn polynomials. We see, that in fact these polynomials coincide with 
the dual -1 Hahn polynomials for odd. Note that in [5], [6] the diagonal coefficients 
vanish 6„ = 0. This corresponds to a simple shift of the argument of the polynomials 
Pn{x) — )■ Pn{x — 1) as is seeu from (14. 2p . 

Consider now the case of N even. The recurrence coefficients are given by formulas 

(ElSD 

The mirror symmetry condition (14. ip is again equivalent to the condition a = /3. 
We then have 

M„ = 4[n]^[N -n + 1]^, 6„ = 4e(-l)"+' - 1, (4.6) 

where 2^ = a — — 1. The positivity condition implies > —1/2- 

The spectrum Xs of the corresponding Jacobi matrix J coincides with the BI grid 
(12. 7p and consists of two uniform grids (with a step of 4 between the neighbor points) 
containing + 1 points and points. There is a gap of length 4 (a — A^) > between 
these sets. 

Condition (i) in this case is equivalent to the restriction 

Ml 

a^7V + ^, (4.7) 

where Mi,M2 are positive coprime odd integers. 

This example seems to be have been overlooked. Note that in contrast to the 
case with A^ odd, the diagonal recurrence coefficients 6„ in (14.61) are not constant. The 
only case when 6„ = — 1 occurs when ^ = which corresponds to the Krawtchouk 
polynomials. 

It is interesting to point out that there is a connection based on the Christoffel 
transform between these two spin models with PST. To see this, suppose we have two 
models A and B. Assume that the spin chain A with PST corresponds to the Jacobi 
matrix Ja with A^+ 1 spectral points Xq, Xi, . . . ,xn that satisfy condition (i). Similarly, 
take spin chain B with PST to correspond to the Jacobi matrix Jb with A^ eigenvalues 
Xq, Xi, X2, . . . , X]si-i also Satisfying (i). The matrix J a has dimension A^ + 1 x A^ + 1 while 
the matrix Jb has dimension N x N and the spectrum of Jb differs from the spectrum of 
J A by the elimination of the level x^. It is then easy to show [llj that the matrix Jb is 
obtained from the matrix Ja by a Christoffel transform. Equivalently, this means that 
the monic orthogonal polynomials Pn{x) corresponding to the matrix Jb are obtained 
from the polynomials Pn{x) corresponding to matrix J a by the formula 

^ ^ P„,.(X)-A-..P„(x)_ 

X — Xn 



where 
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Formula fl4.8p is equivalent to the well known Christoffel transform for the orthogonal 
polynomials [7]. The corresponding recurrence coefficients can be obtained by the 
formulas pT| 

Mn=Mnl7-^, = b^+l + Kn+l - Kn- (4.10) 

Returning to our systems, assume that the Jacobi matrix Ja has size + 1 x + 1 
with odd N and corresponds to the spin chain with coefficients (14. 2p . The spectrum of 
this matrix consists of + 1 points Xg given by (14. 3p . f l4.4p . This matrix corresponds to 
the PST model proposed in The Jacobi matrix Jb has size N x N and its spectrum 
consists of the two subgrids: 

G_ = {Xo, Xi, . . . , X(Ar_i)/2} = G- (4.11) 

and 

^+ = {XiN+l)/2, a^{Ar+3)/2, • • • Xn-i}, (4-12) 

i.e. the subgrid G- coincides with G- and contains (A^ + l)/2 points. The subgrid G+ 
contains (A^ — l)/2 points. 

In order to obtain the recurrence coefficients «„, 6„ corresponding to the matrix Jb 
we need the expression for the coefficients Kn resulting from (14. 9p . A simple calculation 
gives 

_ r2(A^ + a-n) if n even 
" 1 2iN-n) if n odd ' ^ ^ 

Equivalently, in terms of yU-numbers (12. 5p . we have 

Kn = 2[N-n]^/2. (4.14) 
Substituting this expression into formulas (I4.10p we obtain 

Un = 4N,/2[iV - n]„/2, 6„ = -3 - 2(-l)"a. (4.15) 

These recurrence coefficients correspond to (14. 6 p (replacing a — )■ 2^ and performing a 
shift of the coefficient bn)- 

We have thus demonstrated how the spin chain with the PST property for even A^ 
can be obtained from the corresponding spin chain for odd A^. 

It was shown in [11] that provided the eigenvalues Xg, s = 0, 1, . . . of the Jacobi 
matrix J satisfy the condition (i), the weight function Wg constructed from (13. 9 p 
determines uniquely the mirror-symmetric matrix J and hence, the spin chain with 
PST. 

Thus, the only spin chain with the PST property, corresponding to the eigenvalues 
Xg of the BI type (l2.7p .f l2TT2|) are those whose Jacobi matrices J generate the dual -1 
Hahn polynomials with /3 = a. We should stress, that the positivity conditions (12. 6 p 
and (I2.16P are crucial. They imply, in particular that the BI grid consists of two uniform 
subgrids separated by a positive gap. Only for BI grids of such kind can one associate 
the dual -1 Hahn polynomials uniquely. More general BI grids with an overlap of the 
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two subgrids lead to more complicated orthogonal polynomials. This problem will be 
considered elsewhere. 
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